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The role of liquid crystalline mesophases in poly- 
mer crystallization has seen increased attention fol- 
lowing the observations of meta-stable liquid crys- 
talline pre-ordering in the crystallization of isotactic 
polypropylene and other flexible polymers [1-7]. In 
addition to these recent observations, both in Na- 
ture and industry liquid crystal mesophases are used 
to create polymeric materials with desirable proper- 
ties such as spider silk [8] and Kevlar [9]. A bet- 
ter understanding of the theoretical fundamentals be- 
hind formation of the formation kinetics of polymer 
mesophases and their subsequent effect on crystalliza- 
tion would enable many new materials and processing 
techniques to be developed. 

A recent study by Tokita et al [10] on a polymer 
liquid crystal that exhibits a direct isotropic/smectic 
transition has experimentally determined the exis- 
tence of meta-stable nematic orientational order- 
ing preceding the formation of translational smec- 
tic ordering. This was accomplished through the 
study of a polymer liquid crystalline material BB- 
3(1-ME) that exhibits extremely slow liquid crys- 
talline transition dynamics, also developed by Tokita 
et al [11]. The phase-ordering dynamics of this 
material occur at time scales accessible via conven- 
tional methods of polarized light scattering and syn- 
chrotron wide-angle X-ray diffraction analysis tech- 
niques [10,11]. Two key conclusions were determined 
from this work. First, at high-quench rates/super- 
cooling meta-stable nematic (orientational) ordering 



was observed preceding full smectic (orientational 
and translational) order. Second, the occurrence of 
nematic pre-ordering (high-supercooling) resulted in 
morphological changes of growing liquid crystalline 
domains compared to solely smectic growth (low 
super-cooling). Specifically, samples cooled at rates 
high enough to exhibit nematic pre-ordering formed 
well-oriented or "neat" tactoidal smectic domains 
while samples cooled at lower rates, where only smec- 
tic ordering was observed, formed radially oriented 
or textured spherulitic domains. It is noted that tac- 
toids are ubiquitous transient non-spherical drops fre- 
quently observed in phase transformations involving 
mesophases [12,13]. 

The objective of this communication is to show 
that the essential features of the experimental obser- 
vations of Tokita et al [10] described above emerge 
from mesophase phase transformation phenomena 
as described by well-established phenomenological 
Landau-de Gennes type models. 

The existence of meta-stable pre-ordering in phase- 
ordering growth, where the order parameter is non- 
conserved (such as orientation), was first demon- 
strated theoretically by Bechhoefer et al [14], and 
later for dual non-conserved order parameter sys- 
tems [15]. In this past work, Landau-Ginzburg phase 
transition models composed of scalar order parame- 
ters were shown to exhibition a front-splitting insta- 
bility under strong super-cooling and also predicted 
this to occur when large differences in transition time- 




Figure 1: A schematic of a growing liquid crys- 
talline front summarizing the phenomena of inter- 
est: nematic (orientational) and smectic-A (lamellar) 
liquid crystalline ordering and interfacial splitting. 
The orientationally/translationally-order smectic-A 
mesophase is on the left (blue), orientationally- 
ordered nematic mesophase is in the center (green), 
and the isotropic liquid (no orientational or transla- 
tional order) is on the right (red). 

scales exist. 

The scalar models used by Bechhoefer et al [14,15] 
capture the basic pre-ordering phenomena but are 
not suitable for comparison to experimental obser- 
vations. This is especially true for liquid crystalline 
materials where phase-ordering introduces additional 
physics due to inherent anisotropy (see Figure [1]). 
Smectic liquid crystals include dual phase-ordering 
phenomena where orientational order is required for 
translational order to exist. Thus with such a system 
three general effects influence the possible observa- 
tion of meta-stable pre-ordering: (1) the driving po- 
tential (as a function of temperature, concentration, 
and pressure), (2) energetic coupling between the two 
types of order, and (3) phase transition kinetics (dif- 
ference in time scales governed by viscosity ratio). 

In this communication we demonstrate that the 
use of a high-order model and numerical simulation 
is able bridge the gap between the experimental ob- 
servations of Tokita et al [10] and the fundamental 
observations of Bechhoefer et al [14, 15]. As will be 
shown, this approach both captures meta-stable ne- 



matic pre-ordering in addition to actual multi-scale 
dynamic texturing processes. A high-order Landau- 
de Gennes type phenomenological model [16, 17] for 
the direct isotropic/smectic-A liquid crystalline tran- 
sition is used with material parameters based on the 
12CB (dodecyl-cyanobiphenyl) liquid crystal. This 
model does not take into account microstructure that 
is present in polymer liquid crystals, but is an ade- 
quate first-approximation for the system studied by 
Tokita et al [10]. The variation of the time scales 
of nematic (orientational) and smectic-A (transla- 
tional) ordering is easily achieved with this model and 
ideal quenches are assumed. The effects of changing 
the energetic coupling between the orientational and 
translation ordering of the smectic-A mesophase, in 
addition to convection and heat of transition, are left 
to future work in that determination of a suitable set 
of phenomenological constants and model extensions 
are non-trivial tasks [18]. 

Theoretical characterization of orientational and 
translational order is accomplished using order pa- 
rameters that adequately capture the physics in- 
volved. Partial orientational order of the nematic 
phase is characterized using a symmetric traceless 
tensor [16]: 

Q ^ S (^nn ~ ^ij + {mm - II) (1) 

where n/m/1 are the eigenvectors of Q, which char- 
acterize the average molecular orientational axes, and 
S/P are scalars which characterize the extent to 
which the molecules conform to the average orien- 
tational axes [19-21]. Uniaxial order is character- 
ized by S and n, which correspond to the maximum 
eigenvalue (and its corresponding eigenvector) of Q, 
S — ^fin- Biaxial order is characterized by P and 
m/l, which correspond to the lesser eigenvalues and 
eigenvectors, -P = -f {^J■m - M/)- 

The one-dimensional translational order of the 
smectic-A mesophase in addition to the orientational 
order found in nematics is characterized through the 
use of primary (orientational) and secondary (trans- 
lational) order parameters together [22]. A complex 
order parameter can be used to characterize transla- 
tional order [16]: 

^/ = il:e''t' (2) 
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where 4> is the phase and if) is the scalar amplitude 
of the density modulation. The density wave vec- 
tor, which describes the average orientation of the 
smectic-A density modulation, is defined as a = 
V(/)/|V0|. The smectic scalar order parameter ■;/; 
characterizes the magnitude of the density modula- 
tion and is used in a dimensionless form in this work. 
In the smectic-A mesophase the preferred orientation 
of the wave vector is parallel to the average molecular 
orientational axis, n. 

A Landau-de Gennes type model for the first order 
isotropic/smectic-A phase transition is used that was 
initially presented by Mukherjee, Pleiner, and Brand 
[16, 17] and later extended by adding nematic elastic 
terms [23-25]: 



The Landau-Ginzburg time-dependent formulation 
[27] is used to capture the dynamics of the phase 
transition. The general form of the time-dependent 
formulation is as follows [27]: 
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where fir/ fis is the rotational/smectic viscosity, and 
V the volume. A higher order functional derivative 
must be used due to the second-derivative term in 
the free energy equation ([3]): 
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where / is the free energy density, /o is the free en- 
ergy density of the isotropic phase, terms 1-5 are 
the bulk contributions to the free energy, terms 6- 
7 are couplings of nematic and smectic order; both 
the bulk order and coupling of the nematic director 
and smectic density- wave vector, respectively. Terms 
8-10/11-12 are the nematic/smectic elastic contribu- 
tions to the free energy. T is temperature, Tj^j/T^i 
are the hypothetical second order transition temper- 
atures for isotropic/nematic and isotropic/smectic-A 
mesophase transitions (refer to [26] for more detail), 
and the remaining constants are phenomenological 
parameters. 



(6) 

where 9 corresponds to the order parameter. 

As shown below, in this communication we cap- 
ture the essential experimental results [10] by vary- 
ing the quench depth and the translational/rotational 
viscosity ratio. This choice follows from the fact that 
the phase transformation velocity is given by ratio of 
the temperature-controlled free energy driving force 
and the viscosity associated with mesophase order- 
^ ing [28]. This is first studied by treating the kinetic 
phase transformation process as a one-dimensional 
propagation transition front (s), where the unstable 
isotropic phase is replaced by the stable smectic-A 
phase. Multiple one-dimensional simulations were 
conducted for different viscosity ratio/quench depth 
values and the smectic- A/meta-stable nematic front 
behavior characterized. These results were then used 
to determine a dynamic phase diagram in order to 
guide higher-dimensional simulations. Following this, 
a pair of two-dimensional simulations were performed 
in regimes where meta-stable nematic ordering is (1) 
not present and (2) present, guided by the previously 
determined dynamic phase diagram. All computa- 
tional results were achieved by the simulation of eqn 
m using the same computational methods defined in 
ref. [18]. These two stages of the presented work cor- 
respond directly to the two key conclusions of the 
experimental work of Tokita ct al [10]. 
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(a) (b) 



Figure 2: a) Plot of free energy density versus quench depth A = Tai — T for neat smectic-A order (sohd line) 
and degenerate nematic order (dotted line) b) Log-log plot of quench depth versus viscosity ratio = ^ 
(ratio of timescale of orientational and translation order) of the numerical determined SDC (squares/dotted 
line), FDC (circles/solid line) boundaries, and temperature at which the degenerate nematic phase becomes 
unstable (dot-dashed line); see text. The material parameters and phenomenological coefficients, based upon 
12CB [18], arc Tni = 322.85X, Tai = 330.5X, ao = 2 x 10'^;;^, b = 2.823 x 10^ c = 1.972 x 10'';j^, 
ao = 1.903 X lOf';;^^^, = 3.956 x 10^-^, S = 9.792 x 10^ Jj, e = 1.938 x IQ-^^pN, h = 1 x IQ-^^pN, 
6i = 1 X 10-i2piv762 = 3.334 x 10~^°Jm, jin = 8.4 x IQ-^^. 
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Previous work for phenomenological parameter de- 
termination (parameters in eqn[3]) and phase diagram 
computation has been completed for this system [18] 
which shows that a definite free energy minimum ex- 
ists corresponding to smectic-A order. In addition, a 
degenerate nematic-order state exists befow the the- 
oretical second order transition temperature Tai for 
this system. Figure [2^ shows the free energy density 
© versus quench depth (AT = Tai ~ T) for both 
neat smectic-A and degenerate nematic order as pre- 
dicted by the model. As mentioned previously, one- 
dimensional front growth was studied by varying the 
difference in time scales for the formation of nematic 
and smectic-A order, a* = — , and the magnitude of 
the driving variable, temperature quench depth AT. 
Figurel^b shows the computed dynamic growth phase 
diagram where three distinct regimes can be identi- 
fied: 

1. Fully Coupled (FC, lower left region), in this 
regime nematic and smectic-A interfaces cannot 
be distinguished and corresponds to the direct 
transition from an isotropic to smectic-A phase. 
In other words, the nematic pre-ordering layer is 
not observed. 

2. Static Decoupled (SDC, middle region), in this 
regime meta-stable nematic pre-ordering (see 
Figure [2^) is observed to proceed from the ini- 
tial smectic-A front but maintains a constant 
distance/halo and does not fully decouple. In 
this case, a constant thickness nematic pre- 
ordering layer separates the unstable receding 
isotropic phase from the stable advancing smec- 
tic A phase. 

3. Fully Decoupled (FDC, upper right region), in 
this regime meta-stable nematic pre-ordering 
(see Figure is observed to proceed from the 
initial smectic-A front at a velocity greater than 
the preceding smectic-A front velocity. In this 
case, the pre-ordering nematic layer replaces the 
isotropic phase and becomes the matrix in which 
the smectic A phase grows. 

Next we identify the main mechanisms behind the 
results of Figure [2^. The effect of the quench depth 



AT on the kinetics of the phase-ordering transition 
is revealed by the work of Bechhoefer et al [14, 15]. 
Whereas previous work on a dual non-conserved or- 
der parameter system focused on quench depth as 
the driving variable to induce the splitting instabil- 
ity [15] , the effect of differences in dynamic timescales 
of orientational and positional ordering can induce 
the same phenomenon. In the case of low molecular 
mass liquid crystals, as this model is most applica- 
ble, experimental observations point towards the FC 
regime. In the case of liquid crystalline polymers, 
where the time scale for local diffusion increases sub- 
stantially compared to rotation (rigid/semi-rigid and 
side-chain liquid crystalline polymers), the contribu- 
tion of the difference in timescales is predicted to ex- 
hibit this behavior at much lower quench depths, ob- 
served experimentally [10] and in agreement with the 
presented results (Figure [2])b. 

In order to examine the second experimental 
observation, that the presence of meta-stable ne- 
matic pre-ordering induces persistent morphologi- 
cal changes of the liquid crystalline domain, simu- 
lation of the growth of an initially homogeneously 
oriented smectic-A spherulite was performed in two- 
dimensions. This simulation was formed in the FDC 
growth regime (meta-stable nematic pre-ordering) 
and compared to past work [29] on the same system 
in the FC growth regime (no pre-ordering). Quasi- 
bulk conditions were imposed using full periodicity 
which results in an effective ensemble of well-oriented 
spherulites growing in an isotropic matrix. Figures 
[3^ and [2b show the spherulite morphologies at the 
maximum radius achievable with the current compu- 
tational limitations (for the FDC simulation). The 
morphological growth process of the FC regime case 
is complex (refer to [29]), but can be summarized 
as involving the competition between interfacial an- 
choring at the isotropic/smectic-A interface and bulk 
texturing. The resulting spherulite texture assumes 
a quasi-bipolar texture and exhibits an undulation 
instability in the bulk [29]. The effect of the fully 
formed meta-stable nematic matrix phase (FiguretSb) 
has a stabilizing field effect and exhibits tactoidal 
growth of a neat texture. The aspect ratio versus 
long-axis (vertical in both cases) are shown in Figure 
[3t. The effects of nematic pre-ordering are shown 
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(a) (b) (c) 

Figure 3: Two-dimensional simulation results of the growth of an initially homogeneously oriented smectic-A 
spherulite without a) (^* = 25) and with meta-stable nematic pre-order b) {fj,* = 250) where the surface 
corresponds to i?e(^') (the smectic-A layers) and arrows denote the presence of nematic order and the 
director orientation (should be considered headless) c) plot of the aspect ratio versus long axis (vertical) of 
the simulations from a) (circles) and with pre-ordering (triangles). Other simulation parameters were used 
from Figure [2l 



in the current simulation to have a strong morpho- 
logical effect on the smectic-A spherulite, minimiz- 
ing the previously observed interfacial anchoring ef- 
fects and undulation instability. The well-oriented 
nematic matrix phase acts in the same way as a bulk 
electric/magnetic-field [30], which promotes tactoidal 
growth of a well-oriented smectic-A domain, consis- 
tent with the observations of Tokita et al [10]. 

The aspect ratio trend for the FC regime simula- 
tion (Figure 131:) shows a decay, after an initial shape 
dynamic period, of the aspect ratio as the spherulite 
size increases. This decay is explained based upon 
a scaling theory derived from the studying shape 
and director-field transformation of nematic tactoids 
[13,31]. A scaling estimate for the aspect ratio of 
the spherulite was determined to obey the following 
relationship to minimize total free energy [31]: 

- « i^3/5^-3/5-j^-l/5 

r 

where R is the major axis and r is the minor axis, K 
the characteristic Frank elastic constant, r is inter- 
facial tension, and V is the volume of the spherulite. 
Equation ([7]) predicts that the spherulite aspect ratio 



decreases with volume converging to 1, in agreement 
with Figure [31:. 

For the FDC case, the stabilizing nematic pre- 
ordering results in a spherulite with a neat tex- 
ture and thus should obey the Wulf-construction [32] 
which determines surface shape by minimizing the 
sum of total interfacial energy and an ideal undis- 
torted bulk contribution [32]: 

F = / j{r)dA + aV (8) 

J A 

where F is the total free energy of the spherulite, 7 is 
the interfacial tension (a function of position r), and 
a the free energy density of the spherulite bulk. For 
the phenomenological parameters used in these sim- 
ulations, planar anchoring (director parallel to the 
interface) is energetically preferred over homeotropic 
(director perpendicular to the interface). Figures 
[5)d and [31; show that this approximation correctly 
predicts tactoidal growth, maximizing the planar- 
anchored interfacial area. 

The simulated system is one example of the possi- 
ble effects of nematic pre-ordering on the formation 
of a stable smectic-A domain. Past work [33, 34] on 
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the growth of initially homogeneously oriented ne- 
matic spherulites implies another possible scenario. 
In that work, deep quenches were found to result in 
the shedding of four +^ disclination defects in the 
nematic case. This is not morphologically what is 
proposed for the smectic-A case [29]. Thus nematic 
pre-ordering could result in a different smectic-A do- 
main morphology due to difference in defect shedding 
characteristics of the preceding nematic front, which 
would serve as a template. These issues will be in- 
vestigated in future work. 

In conclusion, it has been shown that through 
simulation of a high-order Landau-de Gennes type 
model of the isotropic/smcctic-A liquid crystalline 
transition experimental observations of nematic pre- 
ordering of smectic liquid crystalline transitions [10] 
can be studied. Leveraging this type of high- 
dimensional model has been shown to close the gap 
between experimental observations of this phenom- 
ena [10] and generalized theoretical studies [14,15]. 
Phase transition kinetics results presented show that 
nematic pre-ordering results from both thermody- 
namic potential and dynamic differences in phase- 
ordering timescales (Figure [2)d) , as was observed ex- 
perimentally [10]. Two-dimensional simulation re- 
sults presented explain the experimentally observed 
effect of nematic pre-ordering on growth morpholo- 
gies of smectic domains, where pre-ordering of ne- 
matic domain stabilizes a neat texture in the grow- 
ing smectic-A nucleus, promoting tactoidal growth. 
This work sets the basis for further simulation study 
in two/three-dimensions and extension of the exist- 
ing model [16, 17] to account for the microstruc- 
ture present in polymer liquid crystals. Addition- 
ally, these results, in addition to those of Tokita et 
al [10], determine a possible mechanism for influenc- 
ing the material properties of smectic polymer liquid 
crystals through the induction of meta-stable nematic 
pre-ordering. 
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